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Abstract 

Finding the number 2Hq of directed Hamiltonian cycles in 6-cube 
is problem 43 in Section 7.2.1.1 of Knuth's The Art of Computer Pro- 
gramming ([K nlOj ): various proposed estimates are surveyed below. 
We computed exact value: 

H6=U, 754, 666, 508, 334, 433, 250, 560=61*2^ *217, 199*1, 085, 989*5, 429, 923. 
Also the number Aut^ of those cycles up to automorphisms of 6-cube 
was computed as 147,365,405,634,413,085 
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A Hamiltonian cycle in a graph is a cycle that visits each vertex exactly 
once. Let Hn denote the number of Hamiltonian cycles in n-cube{the graph 
of n-dimensional hypercube). An automorphism of a graph is a permutation 
of its vertex-set preserving its edge-set. Let Autn denote the number of 
Hamiltonian cycles in n-cube up to the group of automorphisms of n-cube. 
Let Weightn denote the number of Hamiltonian cycles in n-cube up to the 
weight. (This another equivalence is well explaned in (citePaOl). 

The number i/„ is given for n < 5 in OEIS (On-Line Encyclopedia of 
Integer Sequences) as the sequence A066037: number of Hamiltonian cycles 
in the binary n-cube, or the number of cyclic n-bit Gray codes ( |S108] ). In 
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fact, Hn = \OHn, where OHn is given for n < 5 by the sequence OEIS 
A003042: number of directed Hamiltonian cycles (or Gray codes) on n-cube. 
Finding OH^ is problem 43 in Section 7.2.1.1 of Knuth's The Art of Computer 
Programming ( [KnlO] ). In Volume 4, exercises, there Knuth also improved 
the general lower bound to Hn > — 0{log'^n)Y" . 

Also, Hn = ^EHn^ where EHn is given for n < 5 by the sequence 
OEIS A091302: number of equivalence classes of Hamiltonian cycles (or Gray 
codes) in the binary n-cube. In fact, Hn is a multiple of y since, representing 
the vertices of n-cube as binary n-sequences, any directed cycle starting from 
the sequence of n zeroes induces a permutation on the n bits, namely the 
order in which they first get set to 1. 

The problem to find Hn was originated by Gilbert in 1958 ( |Gi58j ). 
Perezhogin and Potapov, 2001 ( [PePoOl] ) proved that a„ < Hn < bn for 
n ^ oo, where a„ = e2""'('"W-i+°(i)), 6„ = e2"('"W-i+°(i)). Let M„ de- 
note the number of perfect matchings of n-cube; its values up to n = 6 
are given in OEIS by the sequence A005271. Above bound implies that 
lim^^oo^l^ e [1,2]. Feder and Subi, 2009 r |FeSu09] ) proved that this 

limit is 2; so, it holds Hn = Mn 

For n < 6 the present state of art is given in the Table below. 

For n = 2 and 3 this Table can be easily filled by hand. The values 
Hn, Autn and Weightn were obtained by Parhomenko in 2001 ( |Pa01] ) for 
n = A and by Dejter and Delgado in 2007 ( |DeDe07] ) for n = 5. The value 
Weighty was obtained by Chebiryak and Kroening in 2008 ( |ChKr08] ). The 
exact value Hq = 14, 754, 666, 508, 334, 433, 250, 560, computed by us, as well 
as M| and M|, are given in the Table by upper round. 



n 


Hn 


Mn' 


Autn 


Weghtn 


2 


1 


4 


1 


1 


3 


6 


81 


1 


1 


4 


1,344 


73,984 


9 


4 


5 


906,545,760 


3.471 * 10^1 


237,675 


28 


6 


1.475 * 10^2 


2.667* 10^6 


1,473 * 10i^ 


550 



In order to develop the computation solution, we used the BGL (Boost 
Graph libraries ) the C++ library, developed by Siek, Lie-Quan Lee and 
Lumsdaine ( |SQLOl| ) which involves easy graph construction, implementa- 
tion and the effective parallel computing. The program was written in the 
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Microsoft Visual Studio 2010 developing environment, which includes new 
and effective solutions for developers building applications (both managed 
and native) that take advantage of multiple cores. By the right memory 
handling and process parallelization the computing time lasts about of 6 
months. 

Note that H3=3\, = A\ * (2^) * 7, J/g = 5! * (2^) * 617 * 3,061 and 
He=6\*2^*217, 199*1, 085, 989*5, 429, 923; so, the integer ^ with n=3,4,5,6 
has exactly n-3 odd prime divisors. 

In the next Table we list known upper bounds of Hq. Two last bounds 
are upper rounds obtained from asymptotic bounds ([?] and |FeSu09] ) by 
replacing o(l) by 1 and 0, respectively. The absolute error in the upper 
bound, obtained by Feder and Subi, is going to zero when n — )■ oo and it 
became very small even for the small values of n. So, the absolute upper 
bound as about equal to the exact value of Hn and the computing the Hn 
for n < 7 ,perhaps, is not of big practical importance. 



Known upper bounds of Hq 


upper round 


Dixon and Goodman. 1975 pI)iGo75j 
Douglas, 1977 |Do77j 
Silvermann et al., 1983 |SVS83j 
Clark, 2000 jClOO] 
Perezhojin and Potapov. 2001 [PePoOlj 
Feder and Subi, 2009 |FeSu09j 


1.5 * 10^° 
1.1 * 10^5 
3.7* 10^9 
1.3 * 10^0 
4.1 * 10^4 
2.7 * 10^6 


obtained value of Hq 


1.4* 10^2 



The edge-set of 2n-cube can be partitioned into n Hamiltonian cycles 
( |ABS90j ). To each such partition it correspond 2"^^ Hamilton orientations 
of 2n-cube obtained by orienting one of the cycles and selecting one of 2 
possible orientations on each of remaining n — 1 cycles. In forthcoming paper 
we enumerate Hamilton orientations of 4- and 6-cube and study quasimetrics 
associated with each orientation. 
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